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1.0 INTRODUCTION 

This report summarizes the work accomplished at Rensse- 
laer to study phase distribution phenomenon under simulated 
microgravity conditions. Our group at Rensselaer has been 
able to develop sophisticated analytical models to predict 
phase distribution in two-phase flows under variety of condi- 
tions. These models are based on physics and data obtained 
from carefully controlled experiments that are being con- 
ducted here. These experiments also serve to verify the mod- 
els developed. 


motion of turbulent dispersed solid/fluid flows, the generic 
equation is first multiplied by the phase indicator function, 
Xfc, defined as: 

1, if x e * ^3) 

0, otherwise 

and then averaged over all the possible realizations (ensem- 
bles). Drew et. al. [1995] define ensemble averaging of a func- 
tion F at position x and time t, for some particular realization 
It. as: 


2.0 ANALYSIS 

2.1 Conservation Equations 

The generic conservation equation for each phase is given 

by: 

JLp'p + vp'Pv = v- j + pf (l) 

The corresponding jump condition between components k and 
1 is: 

[(p'F(v-v) + J) ■■]„ - *f, 
where: 

*P ... conserved quantity 
J ... conserved quantity’s flux 
f ... conserved quantity’s source density 
Mj ... conserved quantity’s interfacial source 
Typical state variables are given in Tkble I. 


TABLE I State variables in generic conservation 
equations and jump conditions 
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To obtain the conservation equations which govern the 


F(x,i) = jF(x,iai)<fai(p) (4) 

e 

where E, ensemble, is the set of all the possible realizations. 
Also, if we consider a point on the interface while moving 
with it, we see that the phase indicator function, X k , is a con- 
stant jump discontinuity. Therefore, its material derivative 
must be zero. That is. 


D^JX k 
Dt Tt 


+ T.VJf t «0. 


(5) 


By using the above equation and procedure, ensemble-aver- 
aged conservation equations for adiabatic two-phase flow are 
obtained. These have been described in detail elsewhere [Ala- 
jbegovic 1994], 


2.2 Closure 


2.2.1 The Intarfaclal momentum source 
The interfacial momentum source can be partitioned into 
drag and non-drag components [Lahey & Drew, 1990], as: 

M* = -T VX, = M* (0) (6) 

For drag we may assume: 

M f (D) = -M< 0) = \p e C D \vW (7) 

The non-drag component may be obtained assuming in vis- 
cid flow around spherical particles [Arnold, 1988; Park, 
1992]. The average quantities were obtained using cell-aver- 
aging techniques [Nigmatulin, 1979], which are an approxi- 
mation to ensemble averaging. The resulting expression 
includes virtual mass effects and a lateral force induced by the 
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non-uniform flow field around a spherical particle. The latter 
force is called lift force and is given by [Drew & Lahey 1987; 
1990]: 

M c m c r. e <<PcV xVx, f (8) 


and negligible collisions and small gradients in the dis- 
persed phase's mean velocity, we obtain 





c.j 


(15) 


The efTects of the liquid turbulence interaction with parti- 
cles/bubbles are included in the turbulent dispersion force: 


M <w0).ro 




C TDPc*c^ a d 


(9) 


2.2.2 Binary Collision* botwssn ths particles and the 
particle-wall collision* 

Ahybegovic [1994] has developed closure laws for the 
stress induced in the dispersed phase and the transfer of the 
kinetic energy. For brevity only the procedure used to derive 
the expressions is described here. The following assumptions 
were made: collisions between more than two particles at the 
same instant are negligible; all collisions are elastic with neg- 
ligible rotational effects. Starting from the expression for the 
average stress inside the dispersed phase given by Bachelor, 
[1970] an integral expression fen- the collision-induced shear 
stress is obtained by evaluating the probabilities of collisions 
of particles moving with given velocity, position and size. 
This expression is similar to that obtained by Chapman and 
Cowling [1970]. This was further integrated to obtain the final 
expressions. Similar procedure was used to calculate particle- 
wall collision-induced force. 


2.2.3 Turbulent Interfacial Work 
The interaction between the dispersed particles and the 
eddies is governed by the following equations: 



(v' -v' .) 

x * p, ‘ e '‘ 


( 10 ) 


5(*/ c + VV -° < n > 

The interfacial work due to particle interactions with the 
continuous phase turbulence is: 


The above models have been tested against experimental 
data obtained at Rensselaer and elsewhere for various flow sit- 
uations, and show remarkable accuracy. 

3.0 Solld/Llquld Experiments 

In a two-phase mixture for low volume fractions, the lift 
force is one of the most important phenomena affecting the 
radial distribution of the flow. This effect is due to the relative 
velocity between the two phases. Serizawa et al. [1988] dem- 
onstrated that this lift force pushes bubbles towards the wall in 
upward flows and towards the center of the conduit for down- 
ward flows. 

Figure- 1 shows how the gravity affects the hydrodynamic 
behavior of a turbulent two-phase flow. We may eliminate the 
buoyancy force by using a dispersed phase of same density as 
that of continuous phase. Assad et. al. [1995] measured the 
liquid and dispersed phase distribution in a solid/liquid two- 
phase flow in a horizontal pipe using a Laser Doppler Ane- 
mometer (LDA). The dispersed phase consisted of -2.0 mm 
average diameter size expanded polystyrene particles with 
specific gravity of 1.03. The test section was 30.6 mm in 
diameter and the pipe was made out of a special optically 
clear material (FEP, Fluorinated Ethylene Propylene) with the 
same index of refraction as water. Measurements for the 
velocity profile, turbulent structure and the volume fraction 
for each phase were made for various flow rates and global 
volume fractions. 

The two-fluid model described above was incorporated into 
the PHOENICS code to test it against the data for neutrally 
buoyant particles. Although the model was developed for ver- 
tical flows, only, the neutral buoyant particle data for horizon- 
tal flows can still be compared against the calculations. Figure 
3 shows that the agreement with data is good. 


W" * c W v’ v ' - C*V v' ) 

4. 9 T. V C 1 v «. i V e,j C 2 V pO, > v pO.j) 


where. 


4P„ P, 

* = 3 P'C D \V p -V c 


( 12 ) 


(13) 


2.2.4 Algebraic Stress model for the dispersed phase 
Reynolds Stress 

Starting with the assumption that: 

v'.v'. 

-j* s const. (14) 


4.0 Llquid/Llquld Experiments 

The particles which were used by Assad [1995] were solid 
spheres with a non -deformable, non-slip interface. The flow 
dynamics around them are well understand and documented. 
In case of bubbles and droplets, there is no solid interface; the 
fluid outside interacts with the fluid inside giving rise to inter- 
nal circulations and deformations. Pressure differentials can 
deform the surface so that they are no longer spherical. All 
this affect the interfacia] momentum transfer and the turbu- 
lence structure in a significant way. However, if the droplets/ 
bubbles are small enough, the deformations is very small 
compared to their size and hence can be neglected. The liquid/ 
liquid experiments allow us to further verify our models and 
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collect data to guide us in further enhancements. 

Butyl benzoate was chosen as the dispersed phase for simu- 
lating microgravity. It has a specific gravity of 1.0 at 20° C. 
Also, it is more sensitive to temperature as compared to water. 
This makes it convenient to gravity separate the oil from 
water by heating the mixture by 20°C. It has since been shown 
to have satisfactory optical properties for Phase Doppler mea- 
surements. 

The closed loop mode of operation was found to be unsuit- 
able for obtaining mono-dispersed droplet flows. This moti- 
vated us to consider open loop operation with an injector to 
produce a mono-dispersed suspension. A droplet injector was 
designed and has been constructed with the following charac- 
teristics: (i) ability to control the flow rate of water, (ii) ability 
to control the diameter of the oil droplets, (iii) ability to con- 
trol the relative flow rates of water and oil to obtain desired 
global void fraction, and (iv) most importantly, the droplet 
size distribution should be as small as possible. Indeed, we 
want a mono-dispersed flow so that a correlation between the 
size of the droplet and the phase distribution phenomenon 
may be established. The obvious choice to achieve the above 
objectives is to have capillary tubes as injection units. A pre- 
liminary calculation was done to predict the size of the drop- 
lets as a function of oil volume flow rate, velocity of the 
stream trying to shear off the droplets from the tips and the 
size of the capillary. This was guided by the work done by 
several researchers [Hayworth et al., 1950; Null eL al., 1958] 
where a balance between the drag on the spherical droplet 
being formed at the tip due to the fluid motion and the cohe- 
sion forces due to surface tension yielded a correlation 
between the droplet size form and the capillary size with fluid 
velocity as a parameter. This droplet injector has been tested 
and appears to be functioning as predicted. 

During initial tests, LDA data was found to be inadequate 
to estimate the size of the droplets accurately. This inadequacy 
stems from the method of interpreting the transit time of a 
burst as an indicator of the size of the droplets. This mode of 
droplet size estimation can only be used to get an approximate 
idea of the size distribution. For more accurate data we may 
use a Fiber Phase Doppler Anemometer (FPDA) which relies 
on the phase difference between the reflected light detected by 
two detectors, a fixed distance apart, to determine the size of 
the droplet. A FEP conduit was used as the test section and 
both the probes were submerged under water so that the opti- 
cal path length does not change when traverses are made. 

5.0 The Experimental Facility 

A batch mode oil/water loop was built and tested as shown 
in figure-2. To avoid the realignment of the FPDA probe at 
each measurement location, the probes and the FEP test sec- 
tion are all immersed in water so that no correction or align- 
ment has to be done for measurements at different locations. 

The oil droplet injector was tested to find out the effect of 


oil and water flow rates on the size distribution of droplets. A 
small test loop was constructed outside of the main test facil- 
ity for ease of observation and control. Following conclusions 
were drawn from the data obtained: 

(1) . The droplet size decreases with an increase in water 

flow rate. 

(2) . The average droplet size increases with the increase in 

oil flow rate and at the same time the spread in the size 

distribution increases. 

6.0 Results and Discussion 

A study was done to verify that an FPDA can measure oil 
droplet size. This was necessary to optimize the setup vari- 
ables and verify the optical properties of oil in the context of 
the FPDA measurements. Figure-4 below shows the experi- 
mental facility used for this purpose. A small glass tank was 
filled with a known amount of water. It was determined that 
the reflection mode of diameter measurement is most suited 
for actual measurements. Two void fractions were measured: 
0.07% and 0.21%. Larger void fractions could not be tried 
because of the formation of an emulsion by the stirrer. This 
increased the turbidity of the mixture to such a large extent 
that the measurement rate decreased by a factor of 10. 

Figure-5(a) shows the seeding size distribution in which 3.5 
cc of the seeding was used to seed 19.6 liters of water. A sam- 
ple of the seeding was analyzed under a travelling microscope 
to determine the sizes present. The glass seeding sizes were 
found to be in the range 5 pm to 100 pm, with most of them in 
25 to 50 pm range. Although the data is not statistically repre- 
sentative of the distribution, it gave us the bounds on size. The 
histogram in figure-5(a) is accurate to that extent with most 
particles the same range as above. Optimization of the align- 
ment and different parameters led to different trials with most 
of them resulting in the similar distribution. 

Figure-5(b) shows the particle size distribution in unseeded 
water. The particles measured are the natural seeding sus- 
pended in the water supply. Most of the particles are in range 
0 to 20 pm. Figure-5(c) shows the oil droplet size distribution 
for a volume fraction of 0.07%. Immediately after the mea- 
surements were done with unseeded water, 15.5 cc of oil was 
added. The stirrer immediately broke up the oil globules into 
small droplets and a delay of about 10 minutes was allowed 
before measuring to allow the system to stabilize. The histo- 
gram in figure-5(c) shows die range of sizes present is 0 to 
340 pm. Most of the particles are below 100 pm in diameter. 
If this is compared with the figure-5(b), it is practically impos- 
sible to say which of the data belong to the droplets and which 
to the natural seeding especially in the 0-25 pm range. This is 
the case with low volume fraction, however al higher volume 
fraction, the data percentage in this range must be small to be 
able to neglect the influence of natural seeding. This was con- 
firmed by adding more oil to the tank, the resulting void frac- 
tion being 0.2%. Figure-5(d) shows a typical data set obtained 
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at this volume fraction. Here most of the droplets are in range 
40 to 100 pm and the fraction of particles/droplets less than 20 
Jim is very small. With the injector in the main loop, we 
should be able to keep the global void fraction relatively large 
enough without any emulsification. Thus we are confident that 
with large sample sizes, the data will be representative of the 
true droplet size distribution. 

7.0 Proposed future work 

The main thrust of the future work is to put an experiment 
aboard the space shuttle. The objective of the experiment is to 
measure phase separation in an air/water mixture. The empha- 
sis is on simplicity and efficiency in conducting the experi- 
ment and to provide adequate instrumentation of acquire a 
complete set of data. These data can then be used for further 
understanding the dynamics of bubbly flows in space. The 
instrumentation will consist of: (i) a two-dimensional pen 
Laser Doppler probe; (ii) a Phase Doppler Anemometer, (iii) a 
real time Fast Fourier Transform system, a phase discrimina- 
tor, a high voltage generator and a demodulator, all fitted on a 
single PC card; (iv) a one-dimensional hot film anemometer, 
(v) a traversing mechanism, and (vi) a CCD camera. We will 
measure two-dimensional velocities with the LDA, one 
dimensional velocities and the local void fraction with the hot 
film probe, size distribution of the bubbles with (he CCD cam- 
era and the FPDA. All these will be measured as a function of 
radial position for two axial locations: one close to the inlet 
and one close to the outlet These data will allow us to com- 
pute bubble coalescence and breakup as well as the develop- 
ing phase distribution. These instruments will be 
progressively tested and calibrated in experiments in drop- 
towers and parabolic flights in Lear Jets. The loop itself will 
be a recirculating loop. A sine pump will circulate the mix- 
ture. For monodispersed experiments, we will separate the 
water from air using a centrifuge. For a arm length of 10 cm 
we will need about 600 rpm to simulate normal gravity. A spe- 
cially designed injector will be used to generate the bubbles. 
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Shear-induced turbulence 



Relative velocity 


Interplay of gravity and buoyancy 


Fig. 1. Parametric dependence of two fluid models 



Fig. 2. Schematic of the new facility. 
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Fig. 3. Comparison between measured and predicted velocity and turbulence quantities for neutral particles. 
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Screen 




(a) With glass seeding 


(b) Unseeded Water 



(c) Oil droplets, Global void fraction=0.08% (d) Oil droplets. Global void fraction=OJ 1 % 

Fig. 5. Particle/droplet size distributions. 
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